Fluid Models for Kinetic Effects on Coherent Nonlinear Alfven Waves. II. Numerical 

Solutions 
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The influence of various kinetic efi'ects (e.g. Landau damping, difi'usive and collisional dissipation, and 
finite Larmor radius terms) on the nonlinear evolution of finite amplitude Alfvenic wave trains in a finite-/? 
environment is systematically investigated using a novel, kinetic nonlinear Schrodinger (KNLS) equation. 
The dynamics of Alfven waves is sensitive to the sense of polarization as well as the angle of propagation 
with respect to the ambient magnetic field. Numerical solution for the case with Landau damping reveals 
the formation of dissipative structures, which are quasi-stationary, S-polarized directional (and rotational) 
discontinuities which self-organize from parallel propagating, linearly polarized waves. Parallel propagating 
circularly polarized packets evolve to a few circularly polarized Alfven harmonics on large scales. Stationary 
arc-polarized rotational discontinuities form from obliquely propagating waves. Collisional dissipation, even 
if weak, introduces enhanced wave damping when /3 is very close to unity. Cyclotron motion effects on 
resonant particle interactions introduce cyclotron resonance into the nonlinear Alfven wave dynamics. 
PACS numbers: 96.50.Ci, 52.35.Mw, 52.35.-g, 95.30.Qd 
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I. INTRODUCTION 

The envelope dynamics of nonlinear Alfven waves at 
small-/? are thought to be governed by the derivative 
nonlinear Schrodinger (DNLS) equation, which describes 
parametric coupling with acoustic modes Q. It is be- 
lieved that the dynamics of such waves changes drasti- 
cally in a coUisionless, finite-/^ plasma, due to several 
kinetic effects, especially coUisionless (Landau) damp- 
ing of ion-acoustic oscillations In a highly dissi- 
pative regime, the ion-acoustic gwasi-mode is no longer 
a plasma eigenmode. Thus, it is more natural (accord- 
ing Refs. §,|) to view the mechanism of Landau damp- 
ing of Alfven waves as particle trapping by both the 
magnetic mirror force and electric field. On account 
of the inevitable mathematical difficulties of a kinetic 
approach, many numerical simulations have been per- 
formed. DNLS-bascd simulations with beam pumping 



(comet modeling) 
particle simulations 






hybrid magnetohydrodynamic- 
hybrid particle ions-fluid elec- 
have been performed. 

we systematically accounted 



trons simulations 

In a previous paper 
for ion kinetic effects (but neglected finite conductiv- 
ity, discussed elsewhere [^^) on parallel and slightly 
oblique propagating waves in homogeneous plasma, to 
the order at which the DNLS is derived. A fluid repre- 
sentation for Landau damping allowed us to derive a 
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relatively simple evolution equation for finite amplitude 
Alfven waves in finite-/? environments, typical of the so- 
lar wind. Collisional damping was also included. For 
slightly oblique propagation, the effects of finite Larmor 
radius result in nonlinear dispersion-type integral terms. 
The terms which describe kinetic effects are usually inte- 
gral (nonlocal) in nature, and reflect the finite ion transit 
time through the envelope modulation of an Alfven train. 
However, for the weakly collisional (viscous) dissipation 
and for the small Larmor radius limits, these integral op- 
erators asymptote to (local) differential operators. This 
is not the case for Landau damping, since no intrinsic 
scale exists. The envelope evolution equation obtained 
is thus a nonlocal, integro-differential equation which is 
not not amenable to analytical solution. 

In this paper, we explore numerical solution of the 
kinetically modified derivative nonlinear Schrodinger 
(KNLS) equation. Some preliminary results of this in- 
vestigation are given in Ref. We again classify the 
solutions by the kinetic effects which are taken into ac- 
count. For simplicity, we usually consider parallel prop- 
agating waves, however the results are valid for oblique 
waves since they still obey the KNLS equation ||l^,|l^,|l^ . 

We discover that (especially for the case of Landau 
damping) new types of waveform emerge. The known 
observed "final states" of nonlinear wave steepening (i.e. 
shock structures) are classified as: 

1. ) collisional (hydrodynamic) shocks for which non- 

linear steepening is limited by collisional (viscous) 
dissipation, which sinks energy from small scales 
(i.e. from high-fc harmonics); 

2. ) coUisionless shocks (common in astrophysical 

plasma) in which nonlinear steepening is limited by 
dispersion resulting in soliton-type structures with 
sufficient energy content in high-A: harmonics. 
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A new, large class of waveforms, dissipative structures, 
appears when dissipation acts on all scales (recall Landau 
damping has no intrinsic characteristic scale). Thus, we 
add 

3.) dissipative structures, for which nonlinear steep- 
ening is limited by collisionless (scale invariant) 
damping. Harmonic spectra of such structures con- 
tain predominantly low-fc harmonics, i.e. energy 
resides at small scales. Fast rotation of the of wave 
magnetic field vector (i.e. rotational or directional 
discontinuity formation) is another common feature 
of these structures. 

Dissipative structures emerge in plasma with Ti ~ Tf. for 
a wide range of (3: 0.5 < (3 < 1.5, which is typical for the 
solar wind plasma. Thus, we argue that such structures 
are likely a common constituent of solar wind Alfvenic 
turbulence. Note here that the dissipative structures 
discussed here emerge without external energy source in 
the system. Thus, the amplitude of them decreases with 
time. 

Numerous in situ observations of the solar wind mag- 
netic activity have revealed the nonlinear nature of MHD 
waves 21 . Recent observations indicate the existence 
of directional (e.g. rotational) discontinuities, i.e. re- 
gions of rapid phase jumps where the amplitude also 
varies ||2^,|2^. Several types of directional/rotational dis- 
continuities which might be distinguished by their phase 
portraits have been observed. There are (i) discontinu- 
ities of the "S-type" , at which the magnetic field vector 
rotates first through some angle (less than or close to 
90°) in one direction, followed by rotation in the oppo- 
site direction through an angle larger than 180° (typi- 
180° < < 270°) 



23 1, and (ii) arc-polarized 



cally, 180° < < 270°) [g2| 
discontinuities, where the magnetic field vector rotates 
along an arc through an angle less than 180° |2^,Q. At 
directional discontinuities, the fast phase jump is accom- 
panied by moderate amplitude modulation {5B ^ B). 
At the rotational discontinuities, the amplitude modula- 
tion is small or negligible {SB B). It is shown in this 
paper that all these commonly occurring wave structures 
are manifested in the nonlinear evolution of Alfven waves 
in finite f3, isothermal plasma for different initial condi- 
tions and propagation angle. 

The theory of nondissipative Alfven waves governed 
by the conservative DNLS equation predicts nonlinear 
wave steepening and formation of waveforms with steep 
fronts. These waveforms are unstable to modulations. 
Thus, spiky, many-soliton structures are emitted from 
the steep edge. These oscillatory, (roughly) circularly po- 
larized soliton-like structues interact nonlinearly to yield 
strongly turbulent, spatially irregular wave profiles [see 
next section]. However, the DNLS theory fails to explain 
the formation of rotational and directional discontinu- 
ities. A further step was to include the linear damping 
of Alfven waves due to finite plasma conductivity. It was 
shown that the nonlinear wave relaxes to a shock train 
and constant-B rotational discontinuities that rotate the 



field through exactly 180° ||25|,|26|,|1[ . Inspite of this, the 
DNLS theory was unable to explain the existance and 
dynamics of both (i) the S-polarized directional and ro- 
tational discontinuities and (ii) arc-polatized rotational 
discontinuities with rotation of less than 180°. Recent 
particle code simulations |pT|-jl3| have shown that Landau 
damping may drastically change the dynamics of Alfven 
waves and result in the emergence of such wave struc- 
tures. 

In this paper, we use a recently developed ||l^ ana- 
lytical model of the KNLS equation to investigate the 
influence of Landau damping and several other kinetic 
effects on the dynamics of Alfven waves. The main claim 
of this paper is that all the discontinuous wave struc- 
tures discussed above are distinct solutions of the same 
simple analytical model for different initial conditions, 
e.g. initial wave polarization and wave propagation an- 
gle. A (quasi-) parallel, initially circularly polarized, am- 
plitude modulated wave evolves to a single, purely cir- 
cular, lowest-fc harmonic decoupled from dissipation. A 
(quasi-) parallel linearly polarized wave forms a dissipa- 
tive structure which is characterized by S-type, fast phase 
rotation and moderate amplitude variation, i.e. has the 
requisite properties for classification as a localized, S- 
type directional/rotational discontinuity. Such discon- 
tinuities are intrinsically dissipative, as their amplitude 
continuously decreases with time. An obliquely prop- 
agating wave forms an arc-type rotational discontinuity, 
which does not experience even minimal dissipation. The 
effects of finite Larmor radius (significant at oblique an- 
gles) are shown to tend to suppress arc-polarized struc- 
ture formation. Low-fc, finite-amplitude circular waves 
form instead. CoUisional (diffusive) dissipation of a pon- 
deromotively driven acoustic mode is significant (even 
in a weakly coUisional plasma) if the /3 value is close to 
unity, thus increasing plasma temperature in this region. 

The remainder of this paper is organized as follows. In 
Section ||, we discuss the method of calculation, initial 
conditions and standard DNLS solutions (i.e. the limit- 
ing case of /3 = 0). In Section IH, we solve the KNLS 



for the case of Landau damping. In Section IV, we in- 
vestigate an effect of coUisional damping. In Section ^, 
the influence of finite Larmor radius corrections on the 
dynamics of slightly obliquely propagating waves is in- 
vestigated. Section VI presents a summary of the result 



obtained and some concluding remarks. 



II. NUMERICAL METHOD, INITIAL 
CONDITIONS, AND DNLS SOLUTIONS 

The general KNLS equation may be written as 



db 



1 d 



0, 



(1) 



where b — b{z — VAt,tf,) — {b^ + iby)/BQ is the wave 
magnetic field, va and fli are the Alfven speed and pro- 
ton ion-cyclotron frequency, and tg = {b"^ / B^Y represents 
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slow envelope evolution. The nonlinear ponderomotive 
plasma velocity perturbation is 



ci 



(|6p-(|6p))+C2£[|5| = 



^>]). (2) 



Here Ci and Ci are functions of plasma parameters, i.e. /3, 
Te/Ti, etc., only. Throughout the paper, we set Tg — Ti 
for simplicity. K. is an integral operator which represents 
kinetic effects. The detailed form of it depends on the 
effects under consideration. We now introduce the di- 
mensionless coordinate and time respectively as C, = z/z 
and r — ie/f, where z — bOc/cUp and t ~ 200/51^. Here ujp 
is the proton plasma frequency and c is the speed of light. 
Then, writing b{(,T) = L~^/^ 6^ exp(iAfeC): where 
Xk — 2TTk/L (L is the dimensionless domain length), the 
equation for k-th harmonic, bk, follows as 



ir]\lbk 



dh 
dr 

+ tXk{ci{\b\' 



(3) 



'))b+c,jc[\b\'-{m] b]^ = o, 



where rj = 0.04 is the dispersion parameter. This equa- 
tion was solved for periodic boundary conditions using a 
predictor-corrector scheme and a fast Fourier transform 
to calculate nonlinearities. 

For /3 = no kinetics impacts the wave dynamics [ p^ , 
so Eq. reduces to the familiar DNLS equation with 
Unl simply 



Unl 



1 



VA 



2 1-/3 



(4) 



The DNLS is integrable and has the exact (soliton) solu- 
tion |2^^- |29| jl5[| . We pick a standing (in the wave frame) 
soliton solution (6o is the soliton amplitude): 



e(C) 



bo cxp(i0) 



L/2)bl/2r^y 
arctan {sinh [(C - L/2) bl/2ri] } 



cosh [(C 
2 
3 



(5a) 
(5b) 



To avoid effects of periodic boundary conditions, we have 
chosen the domain length L equal to 327r and taken 8192 
harmonics (and spatial points) so that —4096 < n < 
4096. The test run has shown excellent agreement with 
the analytical solution during the time of computation 
(up to T = 40, i.e. 8,000 cyclotron periods). 

As high-amplitude magnetic perturbations in the so- 
lar wind probably evolve from small-amplitude (linear) 
ones, it is reasonable to examine the nonlinear evolu- 
tion of finite-amplitude periodic waves of both linear and 
circular polarizations. In these cases, we have taken 
L ^ Att and 1024 harmonics and spatial points, i.e. 
—512 < n < 512. Results of reference runs for amplitude 
modulated linear and circular polarizations are shown in 
Figs. |l] and H, respectively. [Note here that only ampli- 
tude modulated waves experience nonlinear evolution, a 
purely circularly polarized wave with b ^ exp(iAC) is an 



exact solution of the DNLS (and KNLS) for which nonlin- 
ear effects do not enter since = 1.] The initial wave 
profiles are given by two initially excited Fourier har- 
monics. For linear polarizations, we pick 6_i — bi — I, 
all others are zeroes, for circular polarizations we pick 
6_2 = fe-i = 1, for elliptical polarizations, we pick 
b-i — 1.1, bi — 0.9. Thus, the waves are left-hand po- 
larized. However amplitude profiles, |6| vs. C, (Figs. ||a, 
||a) look alike, their phase portraits (Figs. ||b) and 
harmonic spectra (Figs. ^e) differ significantly. All 
the waves exhibit the nonlinear steepening phase of a 
front at earlier times (r ~ 2). Dispersion further lim- 
its steepening and produces (at times t 5) small-scale 
parasitic, oscillatory, circularly polarized (even for initial 
linear polarization) wave structures. These oscillations 
are just DNLS solitons radiated from the steep front of 
the wave train by modulation instability. (By different 
authors, such a process is referred to as "dispersive steep- 
ening".) The direction of polarization is shown by an 
arrow on graphs indicating right hand polariza- 

tion. At later times, r 40, nonlinear processes result 
in a wave magnetic field which is completely irregular, 
as seen from Figs. |l|c, |^d, ||d, as well as from their 
harmonic spectra, (Figs. ||c). In these cases strong, 
large-amplitude Alfvenic turbulence is developed. 



III. COLLISIONLESS DISSIPATIVE ALFVEN 
TRAINS 

A. Numerical solutions 

The KNLS equation with coUisionless (Landau) dissi- 
pation [||-|jl|J^,|l[|is Eq. (|l|), together with the nonlin- 
ear velocity perturbation written as 



Unl - 



where 



VA r 

2 I 



Ml 



'-))+M,C[\b\^~{\b\^)]} 



(6) 



Ml 



(l-/3)+xf (1-/3/7) 
(l-/3)2-^xjj (1-/3/7)2' 
(7-l)/7 



^11'^ (1-/3)2-^x^1-/3/7)2' 

7 = 3, c2 = jT/m, 13 = cl/v\, and X|| = (8/3/^7)'/' is 
the parallel heat fiux dissipation coefficient which models 
the intrinsically kinetic effects of resonant particles. Here 
C is the resonant particle integral (Hilbert) operator: 



1 



^/(.')d.' = .A;,. (7) 



The sign = denotes equivalence to a Fourier image and 
V denotes a principal value integral. 
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It is interesting first to investigate the influence of Lan- 
dau damping on a (previously) stable envelope soliton 
wave train. Figs. ||a-^ depict the KNLS waveforms, 
hodographs, spectra at different times as well as wave 
energy degradation and a snapshot of fluctuating fields 
bx and 6y at r = 25 for a soliton initial profile in /3 = 1 
plasma (i.e. Mi = 1.5, M2 = -1.63). The effect of the 
coUisionless dissipation is apparent from the decreasing 
amplitude and increasing width of the waveform (Fig. 
^). Energy "inverse cascades" and accumulates in large 
scales (low-fc harmonics), as shown in Fig. |^c. Note that 
the temporal energy evolution can be separated into two 
phases (see Fig. |^), namely a period of rapid damping 
< r < 5 (so that < t < lO^f^,^^) associated with 
drastic spectral modification, followed by an extended 
quasi-stationary period (of duration of, at least, lO^Q^^) 
during which the energy of the waveform decays very 
slowly. 

The time evolution of an initially linearly polarized 
parallel propagating wave for the same plasma parame- 
ters is shown in Figs. |^. Fig. depicts wave amplitude 
profiles for four times r = 0,1,5,40. Localized quasi- 
stationary structures are seen to form very rapidly, the 
formation time is t/ ~ 2. Again, the amplitude decreases 
and the width of the structures increases with time, due 
to dissipation. Fig. ^ represents the phase portrait of 
the structures at times r = 1,5,40. The dissipative 
structures exhibit an easily distinguishable "S- shaped" 
diagram. It is worth to mention that this characteristic 
shape is well preserved for a wide range of plasma param- 
eters and may be considered a signature of (coUisionless) 
damping (see also Section IV). The harmonic spectrum 
of the dissipative structures is narrow, with most energy 
associated with low-k harmonics (Fig. |^c). For example, 
for r = 40 the harmonic tail is almost negligible after 
\k\ ~ 10. Energy decay is shown in Fig. No distinct 
phases are seen. 

Figs. ^ show the observed localized structures at r = 15 
in more detail. It is seen (Figs. |a and |b) that regions of 
significant field variations are accompanied by fast phase 
rotation through (approximately) tt radians. However, 
in the regions of negligibly varying |5|, linear polariza- 
tion is preserved. This is easily seen from comparison 
of Figs. 1^ and ||d. At the discontinuity (path A-B-C) 
the magnetic field vector completes a rotation through 
TT radians (Fig. p|d). During the subsequent quescent 
region (path C-D), the magnetic field vector resides at 
the "tip" of the left arm indicating pure linear polar- 
ization. At the next discontinuity the vector returns to 
the initial position, similarly completing a tt radian ro- 
tation, as shown by dashed path in Fig. Thus, the 
KNLS dissipative structures have the requisite proper- 
ties to be identified as localized, rotational/ directional 
discontinuities, as recently observed in the solar wind 
1^,^. Note, however, these KNLS rotational discon- 
tinuities are associated with the regions of varying \b\, 
unlike the conventional definition that |6| — const across 
the rotational discontinuity. It may be seen from compar- 



ison with other cases below that there is no sharp differ- 
ence between a rotational and directional discontinuity. 
One may be transformed into another by changing the 
initial wave polarization and propagation angle. Hence, 
we use both words to denote the S-type KNLS discon- 
tinuity. We should also note the remarkable similarity 
of hodographs obtained by solution of the KNLS equa- 
tion and from full numerical plasma simulations 1 11 1. We 
emphasize the quasi-stationary character of the KNLS 
discontinuities, with waveform shape preserved for thou- 
sands of cyclotron times. The structures exhibit narrow, 
localized spectra (Fig. ||3). The inverse spectrum width 
defines the spatial localization of the discontinuity. Also, 
the small-scale parasitic oscillations typical of the fluid 
DNLS are absent (they are probably inhibited by Lan- 
dau damping at earlier stages of evolution). Such KNLS 
discontinuities occur commonly and are not restricted to 
/3's close to unity. These structures are quite evident in 
a wide interval of (3, of approximately 0.5-0.6 to 1.4-1.6. 
The formation time increases from r/ ~ 2 at /3 ~ 1 to 
Tf ~ 15 at the ends of the interval. Thus, the dissipative 
structures still form at smaller M2, however the forma- 
tion time increases when M2 decreases. Wave profiles 
and spectra still look like those in Figs. |^, |5[ The phase 
diagram, however, changes. Its central part expands and 
the "arms" become thicker, so that the whole diagram 
become similar to that shown in Fig. p7| d for the case of 
collisional dissipation. In all cases the phase discontinu- 
ity is localized at the discontinuity of \b\. 

In Figs. 1^ the corresponding waveform dynamics for an 
initially circularly polarized wave with sinusoidal ampli- 
tude modulation is shown for the same plasma parame- 
ters. In contrast to the previous case, circularly polarized 
waves do not evolve into structures with discontinuities 
(smooth, wave- like forms emerge). Rather, they evolve 
in a few r to a single (almost purely) circularly polarized 
harmonic at the lowest k (Figs. |p, |^). A smooth, al- 
most negligible amplitude modulation is imposed on the 
wave. At these regions, weak, arc-polarized phase irreg- 
ularities occur (see Fig. |^,e). Such phase irregularities 
are intermittent structures living just few t, they carry 
little energy (which is mostly in higher-fc harmonics) and, 
thus, disappear quickly. Energy decay (Fig. ^) is sub- 
stantial during the time of spectrum modification, and 
extremely slow at later times. 

Figs. ^ depict the evolution of an initially elliptically 
polarized wave, an intermediate case between purely cir- 
cular and linear polarizations. Plasma parameters are the 
same as in the previous cases. An elliptically polarized 
wave evolves to a wave structure of slowly varying ampli- 
tude (similar to the circular polarization case) as shown 
in Fig. However, there are sudden phase jumps (by 
TT radians) which are localized at regions of varying wave 
amplitude (typical of linear polarizations) . These phase 
discontinuities (which are the semi-circles in Fig. |^b) are 
separated by extended regions of linear polarization dur- 
ing which a tip of the magnetic field vector resides at 
small arms on the graph. Unlike the linear polarization 
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case, the phase discontinuity is not accompanied by the 
wave amphtude discontinuity. A harmonic spectrum is 
again sufficiently narrow and most energy is associated 
with large scales (see Fig. |^c) . Energy dissipation is weak 
in comparison to linear polarizations (Fig. |^d) . 

Obliquely propagating waves are still described by the 
KNLS equation. However, a new wave field which (for- 
mally) contains a perpendicular projection component of 
the ambient magnetic field should be introduced. As- 
suming the ambient field lies in x-z-plane, we write the 
new field as 



b^{b^ + Bo sin 6 + iby)/Bo. 



(8) 



The nonlinear evolution of the linearly and highly ellip- 
tically polarized waves is strongly sensitive to the an- 
gle between the polarization plane and the plane defined 
by the ambient magnetic field vector and the direction 
of wave propagation. This angle is set by initial con- 
ditions. When this angle is small, the oscillating wave 
magnetic field has a longitudinal component along the 
ambient field. Thus, we refer such waves to as longitudi- 
nal. In the opposite case, the wave magnetic field oscil- 
lates (nearly) perpendicularly to the ambient field. Thus, 
such waves are called transverse. Note that classification 
fails for circularly polarized waves, since a polarization 
plane cannot be defied in this case. 

The time evolution of a transverse wave of initial lin- 
ear polarization is shown in Fig. |^, for the same plasma 
parameters and 8 = 45°. Of course, such a value is 
near the boundary of applicability of the KNLS. How- 
ever, we run such an exaggerated case because (i) it takes 
less computer time, as the nonlinear evolution is faster 
and, thus, (ii) it looks more illustrative than those with 
smaller 0. No significant physics is lost, anyway. The 
wave profiles evolve into two KNLS rotational discon- 
tinuities (as in the parallel propagation case) with the 
typical "S-polarization" in a few r (see Figs. and ||b). 
These discontinuities have different group velocities, so 
some time later (r ~ 15) they start to interact with 
each other, as seen in Fig. Since these discontinu- 
ities are characterized by opposite phase rotations, they 
nearly totally annihilate one another duting the interac- 
tion. The final, residual structure has small amplitude 
(b/Bo ~ 0.1), and is roughly arc-polarized, as seen in 
Fig. ^ at r = 30. This final state is extremely weakly 
damped (see Fig. ||f) because of the very small magnetic 
field perturbation amplitude of the waveform (damping 
in KNLS is nonlinear!). Note here that due to redefini- 
tion of the wave field, Eq. (H) , the energy formally reads 



as E = Eu 



(sinO) , where the last term is not as- 



sociated with a wave, but is just a new energy zero-level 
due to the nonzero perpendicular component of the am- 
bient field. The harmonic spectrum evolution is shown 
in Fig. He. 

Figs, g show the evolution of an obliquely propagat- 
ing, amplitude modulated, circularly polarized wave for 
the same conditions. Although some localized struc- 



tures emerge, they are not connected to phase discon- 
tinuities. Nevertheless, they have different group veloc- 
ities and interact (like in the case of a linear polariza- 
tion) to yield a weakly damped, small amplitude, quasi- 
stationary, purely arc-polarized waveform (Fig. |^b) . The 
amplitude of the waveform is sufficiently small {b/Bo ~ 
0.2), but larger than for linear polarization (note that 
|6|/_Bo — const for such a waveform, as seen from Fig. ^). 
Such arc-polarized structures seem very similar to those 
recently observed in the solar wind and identified with 
arc-polarized rotational discontinuities Wc empha- 
size the almost dissipationless character of such discon- 
tinuities (Fig. ^), inspite of their significant harmonic 
energy content, as seen from the broad spectrum in Figs. 

[lo|e. The energy evolution of oblique circularly polar- 
ized waves (unlike all other cases) exhibits three different 
stages (see Fig. ^): (i) the primary spectrum modifica- 
tion when two spatially localized wave structures emerge 
(r < 10), (ii) interaction of these localized wave struc- 
tures (10 < r < 30), and (iii) formation of a residual 
arc-polarized rotational discontinuity (r > 30). The dif- 
ferent stages of spectrum evolution is also shown in Fig. 

A (quasi-) stationary, arc-polarized discontinuity is 
shown in details in Figs. 10, for r = 40. The discontinu- 



ity is associated with minor (almost negligible) amplitude 
modulation ( Fig. p^). As seen from Figs. [lO|a,b,d, the 
discontinuity is a localized structure of typically a disper- 
sion width (~ va/^i). The magnetic field vector makes 
fast clockwise rotation (see the phase diagram. Fig. p^) 
through slightly less than tt radians (path A-B-C). The 
ends A and C are connected by a part of circularly polar- 
ized wave packet (slow counterclockwise rotation in the 
phase diagram, along the perfect arc C-D-A). Circular 
polarization is indicated by smoothly decreasing phase 



outside the discontinuity (Fig. |lO^). Since \b\ ~ const 
across the discontinuity (as well as for a pure circular 
harmonic), it is nearly decoupled from dissipation (note, 
C[const] — 0, and the second [dissipative] term in the 
KNLS vanishes identically). The harmonic spectrum of 
the discontinuity depicted in Fig. |l0| e is quite narrow 
and of low amplitude (remember, the bo harmonic corre- 
sponds to the ambient field, not the wave). Note the re- 
markable similarity of this solution of the KNLS equation 
to the structures detected in the solar wind and observed 
in computer simulations p^ , |l3|| . 

The time evolution of an elliptically polarized, 
obliquely propagating, longitudinal wave (not shown) is 
very similar to the case of parallel propagation until large 
times (i.e. r ~ 60). At later times (t > 60), an arc- 
polarized rotational discontinuity forms. There the mag- 
netic field vector rotates rapidly rotation through an an- 
gle of less than tt radians. 

The nonliear evolution of the transverse waves is com- 
plitely differennt form that of the longitudinal waves, as 
shown in Figs. |ll|. Such waves quickly evolve (i a few t) 
to an arc-polarized rotational discontinuity (Fig. pT|a,b). 
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The magnetic field vector rotates through an angle less 
than 180° at the discontinuity. Energy dissipation is neg- 
ligible in this case (Fig. Od). Harmonic spectra of these 
wveforms are broad (Fig. |ll|c) 

B. The two mode coupling model 



where uj± = 2Mi|6zp(0)p MiE^ . This solution rep- 
resents resonant energy exchange between the two har- 
monics. The typical time scale is of order of the shock 
formation time and is approximately equal to 

^ iM,E,)-\ (13) 



To qualitatively understand the physical processes 
which underly the behavior of dissipative nonlinear 
Alfven waves (of linear and elliptical polarization), we 
construct a simple two-mode coupling model. We as- 
sume the spectrum consists of three low-fc harmonics 
k = 0,±1 (a finite amplitude fc = harmonic corre- 
sponds to oblique propagation, bo = B±/Bq = sin 8. 
This model is based on the results of numerical solutions 
which show that only two (initially excited) A: = ±1 har- 
monics dominate the spectra. Other higher-fc spectrum 
components are of relatively small amplitude and thus 
can be neglected. We also neglect linear dispersion which 
is significant at the steep fronts, i.e. for broad spectra, 
only. Here we emphasize once again that we do not pre- 
tend to give a quantitatively correct answer, but only 
intend to get some insight into the essential physics. 

We write the wave magnetic field in the form b ~ 
+ bo + 6-e~'^ and substitute it into the KNLS 
equation, Eq. (||), together with Eq. (||). The system 
of evolution equations for each harmonic bj is 

b+ + 2i{Mi + iM2) [\boWb+ + b*_) + b+\b-f] =0, (9a) 
bo - 0, (9b) 
6_ - 2z(Mi - iM2) [\bof{b*+ + 5_) + \b+fb^] = 0. (9c) 

We first consider the case of parallelly propagating waves, 
so that bo — sin 8 = 0. The solution of the second equa- 
tion is simply bo — const. Dividing the first equation by 
the third, we obtain the useful equation 



d5+ __ (Ml + zAfs) b*_ 
dbl ~ ~ {Ml ~ iM2) b*^ 



(10) 



While a general solution can be found, we separately con- 
sider (for simplicity) the cases Mi — and M2 = 0. 



1, No dissipation, M2 ~ 

In the dissipationless (fluid DNLS) limit (i.e. (3^0), 
we have from Eq. ([l0| ) 



|6+(t)P + |6_(t)|2 = E^^ const 



(11) 



which is just conservation of energy. The energy is de- 
fined by an initial condition. Using Eq. (|l^), the system 
of Eqs. (|9a|)-(|9q) can be easily integrated to yield 



2. Strong damping, Mi <C M2 

In the case of strong damping (i.e. /? ~ 1), we neglect 
by the term Mi. Then Eq. (|l0| ) yields 



|5+(r)p - \b^{T)\^ = A = const. 



0) 



Unlike the previous case. Landau damping conserves the 
energy difference associated with different Fourier har- 
monics, so that the +k — fc symmetry of the spectrum 
is preserved. Integration of Eqs. (9a)-(p^) with Eq. (|Tl|^) 
results in (A 7^ 0) 

|6+(r)p = |6+(0)p 



(1 + A/|6+(0)|2) exp (8|M2A|t) - 1 

(14a) 



and (A ^ 0) 



\b+{r)\' = 



\b+{or 

l + 8|M2|r' 



(14b) 



Initially elliptical polarizations correspond to 6+ ^ b^, 
thus 



A|= |6+(0)p-|6_(0)P 



> 0. 



Thus, harmonic amplitudes decay exponentially at larger 
times. 



5± cx e-4|M2A|r 
on the characteristic decay time scale of 

rfLay ~ IM2AI-I 



(15) 



(16) 



By contrast, initially linear polarizations with ~ &_ 
(i.e. A ^ 0) decay as a square root of time. 



dH') 
(0') 



6± « l/VlMalr 
with the faster characteristic time-scale 



b±iT)^b±{0) e 



=pzcj-)-r 



(12) 
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3. Quasi-perpendicular propagation 

Now we neglect cubic terms 6±|fo^P in Eqs. (pa|)-(p^) 
as small compared to &±|f'oP — b±. Then these equations 
read 



b+ + 2i{Mi + iM2)\bo\^{h+ + b*_) = 0, 
b- - 2i{Mi - iM2)|6op(&; + 6-) = 0. 



(17a) 
(17b) 



These yield: 



— ib+ + b*_) + 4^(A^l + ^M2)(&+ + b*_)\bo\^ - 0. (18) 

UT 

It is easily seen that harmonic amplitudes decay expo- 
nentially 



5± (X e 



-4|A/2 I singer 



and the characteristic decay time is 



obi 

decay 



iMzsin^ei 



< r. 



lin 

decay 



ell 

decay 



(19) 



(20) 



Thus, given an isotropic initial distribution, wave excita- 
tion in the quasi-parallel direction is most robust while 
waves in the quasi-perpendicular direction are quickly 
damped. 

Figs. |l^ represent the time evolution of two k ~ ±1 
harmonics for the cases shown in Figs. ^, 0, and |^, re- 
spectively. Fig. |l^ corresponds to linear polarization 
(i.e. A = 0). The harmonics coincide all the time and 
decay according to ~ l/V^". Fig. |l2|b is the case of ellip- 
tical polarization (| A| = 0.4). As seen from this figure, A 
is approximately constant and the harmonics decay expo- 
nentially, contrary to the case of linear polarization. Fig. 
[T^ shows the case of obliquely (9 — 45°) propagating, 
circularly polarized waves. 

Amplitude modulated, parallel propagating circularly 
polarized waves with initially excited harmonics —1 and 
—2 can be considered analogously. However, one should 
consider a spectrum which includes —2 < k < +2 com- 
ponents. This complicates the analytical treatment of 
the problem. The characteristic time, however, can be 



estimated from the elliptical polarization case, Eq. (16) 



Now A ~ |&-2p, and the second (i.e. 6-2) harmonic 



dissipates in a time 



{mib- 



(21) 



After this time, only one 5_ harmonic survives. Thus, no 
significant damping occures. 



C. Discussion 

The sharp contrast between the nonlinear evolution 
of initially linearly (elliptically) and circularly polarized 
parallel propagating wave trains is an immediate con- 
sequence of the unique harmonic scaling of coUisionless 



(Landau) dissipation in the KNLS equation. Integration 
of Eq. (0) with Unl defined by Eq. (^) over all space 
(assuming that b vanishes at infinity) yields (in dimen- 
sionless variables) 



dE 
97 



-Ah 
d_ 

M2 
2 



\b\'§-^cm dc 



2tt 

-00 



k" 



^),„dfc")dc 



(22) 



where E = l^pd^ is the total wave energy, and we 
assumed = 0.(This result have been reported also 

in Refs. |^,Q.) Thus, the rate of energy dissipation is 
controlled by the coupling coefficient, M2 < 0. Also, we 
conclude that higher-fc harmonics are strongly damped 
which is typical for a phase-mixing process (i.e. smaller 
scales mix faster). It is crucial, however, to understand 
that coUisionless damping enters at all k, in contrast to 
hydrodynamic systems where diffusion (viscosity) yields 
dissipation only at large k (i.e. small scales, or steep 
gradient regions). Hence, competition of nonlinear har- 
monic generation (i.e. front steepening) and collision- 
less damping produces a quasi-stationary energy spec- 
trum which is a rapidly decreasing function of fc, with the 
long-time asymptotic behavior dominated by a few low-fc 
harmonics. Since linear polarizations have (initial) spec- 
tra symmetric upon fc — > — fc, they couple more strongly 
to dissipation than circular polarizations do. Moreover, 
interaction of finite amplitude spectrum components nec- 
essarily leads to higher-fc harmonic generation and wave 
steepening. Such interactions are excluded for pure circu- 
lar polarization. Therefore, the circularly polarized wave 
evolves to a single harmonic final state, which is, itself, 
a stationary (and exact) solution of the KNLS equation 
(i.e. it experiences no steepening and minimal damping). 

Since Landau damping enters symmetrically for +k 
and —k spectrum components, it does not change the 
initial symmetry of a spectrum, so that the sense of 
an initial polarization (set by initial spectrum symme- 
try) is preserved (neglecting feeble high- A:, low-amplitude 
asymmetry induced by steepening). Thus, KNLS evo- 
lution preserves initial wave helicities forming direc- 
tional/rotational discontinuities at the edges of the re- 
gions of linear (elliptical) polarization. These direc- 
tional/rotational discontinuities are localized at regions 
of strong amplitude variations, unlike the conventional 
view [01 where the shock and rotational discontinuity 
are well separated and have different group velocities. 
The KNLS directional/rotational discontinuities display 
a well distinguished S-shaped phase portrait which can 
be considered a footprint of the occurrence of strong Lan- 
dau damping in satellite observations. Since the KNLS 
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equation preserves (initial) wave helicity |^,^ , rota- 
tional discontinuities do not form from circularly polar- 
ized waves. Given that circular polarizations decay much 
more slowly than linear polarizations do, the results sug- 
gest that large-scale circularly polarized waves should 
occur more frequently than small scale discontinuities 
which evolve from linear polarizations (given comparable 
initial populations) , especially at large distances from the 
Sun. 

Oblique and quasi-perpendicular, linearly polarized 
waves first evolve to two KNLS directional/rotational dis- 
continuities of opposite phase rotation propagating with 
different group velocities, so that they finally merge and 
annihilate to yield weakly-nonlinear, small amplitude, 
arc-polarized waves. Despite the fact that oblique, circu- 
larly polarized waves do not form KNLS rotational dis- 
continuities, their final state is similar to (but of larger 
amplitude than) that of a linear polarization. Since, 
oblique and quasi-perpendicular (both linear and circu- 
lar polarizations) waves are coupled to the ambient mag- 
netic field B± = Bq sin0, which is usually (much) larger 
than the fiuctuating wave field b, the process of energy 
dissipation is B\/b^ times faster for such waves. Wave- 
forms again evolve to the lowest-A: state which is now the 
single harmonic 60 = B±/Bo (which corresponds to the 
ambient field component, not the wave perturbed field) 
indicating dissipation of all wave energy. Thus, the fi- 
nal state of nonlinear evolution for "off-axis" propagat- 
ing waves is a (very) small amplitude, weakly-nonlinear, 
weakly-damped, arc-polarized Alfven wave packet. Pure 
arc-polarized waveforms (which develop from circular po- 
larizations) have much in common with the structures 
recently observed in the solar wind and identified with 
arc-polarized rotational discontinuities [p4| . Such struc- 
tures experience minimal damping. The comparably fast 
damping of oblique and quasi-perpendicular waves sug- 
gests that most magnetic fiuctuations should occur at 
small angles to the ambient field (given initial isotropic 
fluctuations) resulting in large turbulence levels in par- 
allel and quasi-parallel directions and small turbulence 
levels in quasi-perpendicular directions where Alfvenic 
turbulence is sufficiently suppressed. 



IV. COLLISIONAL DISSIPATIVE ALFVEN 
TRAINS 

A. Numerical solutions 

The evolution equation for coherent Alfven waves is 
again Eq. (|l|) with the nonlinear velocity perturbation 
term written as 

Unl = ^ {Qi - im) + Q2JL - (|6P>] } , 



Oi 



1 



(23) 



P 7-1 



1 



Xl 7 (1-/3/7)2 



7 = 3, and Xjj is the parallel heat conduction coefficient 
which represents diffusive (coUisional) energy dissipation. 
Here Jl is the collisional dissipation (integral) operator 



Mf]{x)^e-^"^ / e-'/i/(x')dx' 



1 



1/L-ik' 



(24) 



where 



where L — (x|j/w^)(l — — /?). For future use we 

introduce k, — z/ the dimensionless dissipation scale, 
i.e. the ratio of a typical wave length to the dissipation 
length. For (3 close to unity we may write approximately 

As in the case of coUisionless dissipation, we first in- 
vestigate the temporal evolution of an (initially) soliton 
wave envelope, Eqs. ([5a|)-(5b), since it is stable in ab- 
sence of any dissipation. Solutions at different times are 
shown in Figs. |l| with k = 100 and (3 = 0.99. We shaU 
see that collisional dissipation is strongest when k ^ 
i.e. when the dissipation scale is comparable with a typ- 
ical wave-length of a train. We again (as with Landau 
damping) see rapid primary spectrum and spatial profile 
modification (Figs. [l3[:,a, respectively) and formation of 
a quasi-stationary solution with negligible energy degra- 
dation (Fig. p^). The hodographs are shown at three 
times r = 0, 4, 25 in Fig. |l3|b. The steady-state profiles 
of bx and by are shown separately in Fig. ^3|e. We should 
mention that the effect of collisional damping is signifi- 
cant only for (3 very close to unity. 

The nonlinear evolution of an initially linear polariza- 
tion (in (3 — 0.99 plasma) is shown in Figs. 14, ij, |l^ for 
weakly collisional (k = 200), intermediate (k — 20), and 
strongly collisional (k — 2) regimes, respectively. Energy 
dissipation is represented in Figs. |l8|. In a weakly colli- 
sional regime (Figs. |lj), a wave experiences steepening of 
the back front, as a result of nonlinear dispersion set by 
the sign of k (or L) which, in turn, depends on plasma (3. 
Since the nonlinear dispersion is positive in this case, it 
competes with the linear (negative) dispersion. The non- 
linear dispersion is most important in regions of a high 
wave magnetic field (since it is proportional to \b ■^), while 
conventional linear dispersion [last term in Eq. (l|)] is im- 
portant where |6| is small. Parasitic small-scale circularly 
polarized oscillations appear on the steep fronts when 
wave steepening is limited by dispersion (Figs. p^a,b). 
These parasitic harmonics, however, disappear quickly 
as a consequence of the fact that high-Zc components are 
strongly coupled to dissipation. The magnetic field vec- 
tor rotates rapidly at the steepened edge of the wave. 
After the phase of nonlinear steepening, a propagating 
shock wave separates from the (standing) phase rotation 
region (i.e. the rotational discontinuity). The width of 
the shock front is controlled by collision length L and 
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its amplitude gradually decreases due to damping (Fig. 
[l^). Such a process is remarkably similar to that first 
discussed by Cohen and Kulsrud [Q. However, it is just 
one limiting case {k oo) of the KNLS wave evolu- 
tion. The spatial width of the rotational discontinuities 
is controlled by dispersion. Hodographs of rotational dis- 
continuities presented in Fig. [l4| d are different from the 
Cohen-Kulsrud case. Here the magnetic field vector b 
makes a counterclockwise rotation through some angle 
(typically 50° — 90°) and than makes a clockwise rotation 
to complete a total 180° phase jump across the disconti- 
nuity. As shown in Fig. harmonic spectra are broad, 
as collisional damping enters at small scales k > k. 

Figs. |l^ correspond to intermediate collisional damp- 
ing (k — 20). Parasitic small-scale, small amplitude os- 
cillations are still emitted from the seepened front (Fig. 
|l5|a,b). Propagating shocks (Fig. |l5|c) are sufficiently 
broad and disappear at a higher rate than those in the 
weakly collisional case. Rotational discontinuities are lo- 
calized at the regions of significant amplitude |6| vari- 
ations, and the typical hodograph takes on an S-shape 
(however, with "inflated", almost semicircular central 
part). The harmonic spectrum is narrower than in the 
K = 200 case (Fig. p^e). The rotational discontinuities 
move slowly to the right (as with Landau dissipation). 
Neither the sign of k, nor the direction of shock wave 
propagation affect the propagation direction of these 
phase discontinuities. More details of this case are shown 
in Figs. |l^ at T = 15. A remarkable similarity of these 
figures to those of Figs. |5| (Landau damping case) indi- 
cates a lack of a characteristic dissipation scale in the 
system, i.e. the dissipation scale is comparable with a 
typical wave-length, k ~ fc. 

The strong collisional regime is represented in Figs. |l^. 
This case exactly corresponds to that with Landau damp- 
ing, since collisional and wave lengths are equal (note, the 
fc-scaling of damping is different). No small-scale har- 
monics emerge, and shock wave generation is suppressed 
by damping (Fig. |l6|a), as typical of Landau dissipation 
S-shaped hodograph displays (Fig. p^b). Most of the 
energy is located in a few low-fc harmonics (Fig. [l6|c). 
For a more strongly collisional regime (k = 0.2, results 
are not shown) the effect of dissipation becomes weaker, 
and wave evolution becomes similar to that of the weakly 
damped wave case. This effect is associated with the 
weakness of the resonant coupling (defined by Q2 coeffi- 
cient) of Alfven and ion-acoustic waves. 

The temporal evolution of wave energy is shown in Fig. 
[l8| a for several collisional regimes, from weak collisional- 
ity (k = 200) to a collision dominated plasma (k — 0.2). 
Fig. [l^ b depicts total wave energy content as a function 
of K [i.e. heat conduction coefficient as (xy)"^] at dif- 
ferent times T = 5, 10, 20. The energy dissipation rate 
increases with coUisionality (i.e. with decrease of k) un- 
til k ~ 1 and decreases for k < 1. The last effect is due 
to the weak coupling of Alfven and ion-acoustic waves 
(the coefficient Q2), which occures when the ion-acoustic 



frequency becomes smaller than its collisional damping 
rate. 

Since the sign of k (or L) defines the sign of nonlin- 
ear dispersion, we also considered the case of /? > 1. 
Traces of typical evolution are depicted in Fig. 19 for 
K — —20 and f3 = 1.01. In this case both linear and non- 
linear dispersions are negative (and do not compete), so 
many small scale oscillations appear at the (sharp) front 
of a wave train (Fig. |l|a,b). This " noisy" harmonic 
generation is possibly related to the modulation instabil- 
ity of steep fronts against small-scale modulations [Q. 
The instability growth rate in the coUisionless regime 
decreases with increasing wave amplitude. In the colli- 
sional regime, the growth rate is amplitude independent. 
Thus, the large amplitude suppression which is the case 
for coUisionless damping is absent for diffusively dissi- 
pative Alfven waves. Further evolution of Alfven trains 
generally follows the previous cases where f3 < 1. 

Circular polarizations evolve similarly to the Landau 
damping case. Despite the emergence of a shock wave 
(especially in a weakly collisional plasma), they do not 
evolve to rotational discontinuities. Instead, they gener- 
ate smooth waveforms with irregular phase rotation. The 
sign of K again defines the sign of dispersion. Its value 
controls spatial scales. The result of dissipative evolu- 
tion is a dominant (lowest-fc), circular harmonic and a 
low-amplitude harmonic background tail, spreading up 
to scales k ^ n. For illustration. Figs. ^ show evolu- 
tion of a circularly polarized wave in a weakly dissipative 
plasma (more typical for astrophysical situations) with 
/3 = 0.99 and n = 200. 



B. Discussion 

A significant difference between collisional (diffusive) 
and coUisionless (Landau) dissipation is the existence of 
a characteristic collisional scale for the former case. As 
usual, it represents a "boundary" between the weakly and 
the strongly damped (i.e. low and hirii-/c, respectively) 
parts of wave spectra. Integrating Eq. (|l|) with Unl from 
Eq. ( p3| ) over space (see Section p^) , we obtain: 

If = -0. {jl m + m'i^j^ m) « 



Q2 



)k\ i-i 



ik 



ik 



dk 



dk 



(25) 



Thus, energy dissipation is controlled by both the coeffi- 
cient Q2 < (which, in turn, controls the coupling of the 
Alfvenic mode to the damped ion-acoustic mode), and 
the dissipation scale k oc (1 — /3)/X|j- Unlike the coUi- 
sionless case, diffusion enters at particular scales k > k. 
As a result, shock waves may form and their steepen- 
ing will be limited by the collisional scale k. The sign 
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of K defines the sign of nonlinear dispersion (e.g. steep- 
ening direction, etc.) and does not affect the dissipa- 
tion rate. Since most spectral energy is associated with 
sufficiently small-A: components, the integral in Eq. (|25| ) 
increases with decreasing k until n is about unity. For 
smaller n the integral is approximately constant. The 
coupling coefficient Q2 decreases with decreasing k (for 
constant (3). The higher the plasma coUisionality, the 
weaker the resonance of Alfven and ion-acoustic waves, 
and thus the smaller Q2 is. Thus, the energy dissipation 
rate is maximum when the dissipation length is compara- 
ble with a typical wave-length. For a strongly coUisional 
regime, ion-acoustic wave damping rate becomes larger 
than its frequency, so that ion-acoustic perturbations are 
no longer a plasma eigenmode. Resonant coupling and 
energy transfer from an Alfvenic mode to ion-acoustic os- 
cillations is weak, despite strong dissipation of the latter. 

One may mention the interesting case of strong damp- 
ing {k ~ ^wave)- III this case, the system evolves as if 
subject to collisionless dissipation, i.e. S-polarized KNLS 
dissipative structures emerge from linear polarizations 
and circular polarizations quickly evolve to a single har- 
monic state. However, such a scale-free regime may occur 
only in strongly coUisional plasma, x| ^ ''^^/^ij with /3 
very close to unity (0.95 < (3 < 1.05). Inspite of its rare 
occurrence in nature, this case allows us to understand 
physical the reasons underlying the uniqueness of colli- 
sionless (Landau) dissipation for Alfven wave dynamics. 

The weakly coUisional regime {k —^ 00) seems more 
typical for astrophysical and space plasmas. Nonlinear 
steepening of a wave front is first limited by nonlinear dis- 
persion, giving rise to radiation of small-scale, circularly 
polarized, noisy parasitic oscillations by the steep front 
(probably due to modulational instability). These oscil- 
lations damp quickly. At later times, the front separates 
into a shock wave moving in the direction of primary wave 
steepening and a waveform the type of which (e.g. direc- 
tional or rotational discontinuity, small-scale phase mod- 
ulations, etc.) depends on the initial wave polarization. 
For initially linear polarizations, this waveform is a rota- 
tional discontinuity (with a weak amplitude modulation 
imposed on it). For circular polarizations, it is an inter- 
mittent, wavy structure (the amplitude of perturbations 
of |6| decreases with time), carrying weak, roughly arc- 
polarized phase irregularities. The rotational discontinu- 
ity (for P < 1) consists of partially counterclockwise rota- 
tion through some angle (typically 50° — 90°) followed by 
clockwise rotation through an angle larger than 180° to 
complete a 180° phase difference across the discontinuity. 
These waveforms have a characteristic spatial scale set by 
dispersion. The propagating shock wave front is much 
narrower and controlled by collisions. The amplitude of 
a shock decreases with time due to coUisional (diffusive) 
dissipation. The residual waveforms are weakly coupled 
to nonlinear dissipation because |6| — const across them. 
In general, linear polarizations are more strongly damped 
than circular ones. 



Given that both linear and circular polarizations are 
initially excited, the weakly coUisional quasi-stationary 
turbulence will consist of smooth waveforms associated 
with rotational discontinuities and purely circularly po- 
larized, high amplitude waves. 

Consider again Eq. (|25| ) in a weak coUisional regime, 
i.e. xjj 0. Then, as k cx (1 — /3)/X|j, k is extremely 
large (and dissipation is very weak) for all values of /3, 
except /3 = 1. The coefficient Q2 oc (3/xl is not sin- 
gular at /3 = 1. Thus, for |/3 — 1| 9^ the dissipation 
rate E ~ Q2/k'^ ^ ^l/^'^ ^ ^^'^ ^ indicating weak 
wave damping. However, for [3 close to unity, we infer 
from Eq. (25) that k ^ and E ^ Q2 ^ l/xy which 
is larger, for weaker coUisionality In space plasma, (3 is 
a typically decreasing function of the distance from the 
sun [i.e. P = f3{R)], starting from P > 1 near the sun 
and decreasing to /3 < 1 further from it. Moreover, this 
plasma is only weakly coUisional. According to the dis- 
cussion above, nonlinear waves propagating out of the 
sun are extremely strongly damped in a narrow region at 
the distance Rd where P{Rd) — 1. This region should be 
easily identifiable by a (much) higher (electron) temper- 
ature compared to elsewhere. The width of the region 
of enhanced dissipation may be estimated from the con- 



dition 



1. Thus, the 



75(l-/3(i?))/x^ 
smaller the coUisionality, the thinner the high tempera- 
ture region. Because of strong damping at (3 = 1, waves 
existing at distances either R > Rd or R < Rd never 
transit through the point R = Rd which is thus can act 
as a "sharp boundary" between these two regions. Thus, 
the level of magnetic field perturbations will be low be- 
yond the point Rd (assuming that the Alfvenic fluctia- 
tions have been generated near the sun), provided the 
width of the dissipation region is large enough. 



V. GYRO-KINETIC (FINITE LARMOR RADIUS) 
EFFECTS ON SLIGHTLY OBLIQUE ALFVEN 
TRAINS 

A. Numerical solutions 

The evolution equation in this case is again Eq. (|^). 
The nonlinear velocity perturbation reads 



NL 



Ah 



(^Ja — J-A 



where 



A^ 



2A 



A2 ( 



A^i + N2C 




(26) 



(l-/3) + [x|/7](l-/3/7) 

(i-py+xiii-Phr ' 
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^ (l-2/?)+x^(l-2/3/7) 



N2 = -2t]Ix\\(3- 



7 a-pr+xi{i~Ph) 



•2 — 



77|l = (pj/2)sin 9, and pi is the ion gyro-radius. [Eq. 
(49) of Ref. [Q contains a typographical error, in that 
the two terms in round brackets were omitted] . Note here 
that can take on both positive and negative values. 
The operator C represents the effects of resonant parti- 
cles. The operator which contains Larmor averaging is 
a composition of two coUisional operators J7l . It can be 
written as 



^(jA-J-A)=i 



(27) 



where we redefine k — z/A. Thus it can assume either 
real or imaginary values. 

In the limits kpi <C 1 and 8^1 (i.e. k/n ^ 1), 
we may expand denominators in the right hand side and 
plug into Eq. (|2^). Then, to lowest order, we write 



Unl 



A^O 




+ [(A^Mi - A^i) 



(A^Ms - N2) Z 



(28) 



First, one may notice that there is no sharp transition 
from the A^ > regime to the A^ < regime. Second, 
gyro-averaging the wave fields results in nonlinear dis- 
persion of waves. It does not introduce additional dissi- 
pation. It may, however, generate (due to modulational 
instability) a large number of strongly damped high-fc 
harmonics. Thus, it may affect the wave damping pro- 
cess. 

Figs. ^ and ^ depict the evolution of linearly and 
circularly polarized initial conditions (respectively) in 
(3 < 1 environment for k — 3000 which corresponds to 
kpi ~ 0.08. The angle of wave propagation is 8 = 45° 
in both cases. There is no significant difference observed 
in the evolution of an oblique linear polarization for zero 
and finite Larmor radii (compare Figs. |^ and respec- 
tively). Quasi-stationary KNLS dissipative structures 
(i.e. rotational discontinuities with |6p ^ const) form 
during the evolution. The width of each is controlled 
by both linear and nonlinear (due to gyro-kinetic effects) 
dispersion. The harmonic spectrum, however, is more ir- 
regular than in the case of zero Larmor radius (Fig. |2l|c) . 
In contrast, an obliquely propagating, circularly polar- 
ized wave behaves differently, and its evolution is closer 
to the parallel propagating case. No arc-polarized struc- 
ture forms from an initially circular polarization. Only 
a waveform with dominantly low-fc circularly polarized 
harmonic forms instead (Fig. p^). 

Finite Larmor radius effects enter the energy dissipa- 
tion rate, as shown in Figs. In this case (3 — 1.3, k = 



50 (i.e. kpi ~ 0.5), and 77 = 0.2. Again, two KNLS ro- 
tational discontinuities of opposite senses of polarization 
form. They interact with each other quickly due to their 
increased width (set by dispersion) , despite the fact that 
their group velocities are similar. The combined effects 
of Landau damping and nonlinear Larmor dispersion re- 
sult in rapid wave energy decay. This is shown in Figs. 
||c,d. 

We should mention another interesting case (results 
are not shown). For (3 < 1 and comparably small |«;| 
(< 50 — 100), solutions of the KNLS equation are un- 
stable, i.e. the amplitudes of several harmonics increase 
rapidly with time. Our program does not allow us to 
follow the full dynamics of such solutions, because a fast 
Fourier transform (used in the program) works correctly 
only when harmonic amplitude vanishes quickly with in- 
creasing k. Note however that the perturbative deriva- 
tion of the KNLS fails in this case, too. The physical 
understanding of this instability is discussed in the next 
subsection. 



B. Discussion 

It is easily seen that finite Larmor radius corrections 
alone do not introduce a damping process of any kind. 
Nevertheless, they can change the energy dissipation rate 
via coUisionless damping Indeed, it can even change the 
sign of the dissipation, as will be discussed later. Using 
our previous results, namely Eqs. ( |2^ ) and (]2^), we may 
now easily write 



dE 
97 



,x ^2 



fc2 



§ Mfcl dfc. 



(29) 



It is very important to emphasize that is not necessar- 
ily positive, but defined by the value of fi (i.e. k, oc 1/A). 
Thus, the last equation may be rewritten as 



dE 



■ Mo 



2^ ,2 

)k\ 



No 



+ A;2 A2 k2 -I- fc2 



k dfc. 



> 0; 



= < 



Mo 



No 



|k2| -fc2 |A2| |k2| -fc2 



Ifcdfc, k2 < 0. 



Thus, nonlinear Larmor corrections for high-/? plasma 
(i.e. k2 > may result in enhanced dissipation, only. 
At low-/? (i.e. /3 < 1,k2 < 0) they may also change 
the sign of dE/dr at least for |fc| > |«;|, i.e. result in 
instability (i.e. pumping of a wave by plasma resonant 
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particles). In the case of (3 less than unity, we replace 
K ^ —i\K\ in Eq. {^t\). It becomes 

2t|^(j^.|A|-J-.|A|) = ^P^ + ^^) (30) 

which reveals singularities of the finite Larmor operator 
at ±fc — \k\. The terms on the right hand side change 
their signs at k = ^\k\. Thus, the sign of Landau damp- 
ing (oc C) changes sign for right hand polarized harmon- 
ics with k < — \k\ (the first term) and left hand polarized 
waves with k > \k\ (the second term). Most small-A: har- 
monics are damped on bulk plasma resonant particles. 
Since plasma particles rotate over their Larmor circles (in 
the perpendicular plane) along with longitudinal thermal 
motion, some large-/c harmonics may experience growth 
due to cyclotron resonance with the particles from the 
tail of the distribution function. Thus, finite Larmor ra- 
dius effects together with the effect of resonant particles 
in /? < 1 environment can drive an ion-cyclotron insta- 
bility (for growing harmonics k > \k\) as well as ion- 
cyclotron damping (fc < |k|), enhanced dissipation (due 
to nonlinear dispersion) in /3 > 1 plasma. Of course, 
the ion-cyclotron instability saturates with time by, for 
example, particle distribution function modification, an 
effect which is not included in this simple model. For 
some values of plasma parameters, however, conventional 
Landau damping may overcome generation of harmonics 
and suppress the instability. 

VI. CONCLUSIONS 

In this paper, the influence of various kinetic effects 
(e.g. Landau and diffusive dissipation as well as effects 
of the gyro-averaging of fluctuating fields on a Larmor 
radius scale) were investigated in detail. A tractable an- 
alytic evolution equation, the KNLS equation, was nu- 
merically solved to study nonlinear dynamics of finite- 
amplitude coherent Alfven waves in a finite-/?, Ti ~ Tg 
plasma, natural to the solar wind. Previous DNLS stud- 
ies failed to describe the evolution of the wide class of 
discontinuous waveforms which were recently observed in 
the solar wind. These waveforms are the (i) arc-polarized 
rotational discontinuities with phase difference less and 
equal to tt and (ii) S-polarized rotational and directional 
discontinuities. Current study shows that all these dis- 
continuous wave structures discussed are distinct solu- 
tions of the same simple analytical KNLS model for dif- 
ferent initial conditions, e.g. initial wave polarization and 
wave propagation angle. The results are summarized in 
Table|. 

I. Collisionless (Landau) dampmg introduces a unique 
scale invariant dissipation process into the dynamics of 
Alfvenic trains. In contrast to viscous/diffusive dissipa- 
tion. Landau damping is not restricted to small scales. 
Thus, it is not surprisingly that the harmonic content of 
collisionless dissipative waveforms tends to be dominated 



by large scales, with most wave energy contained in the 
low-fc harmonics. Since there is no characterictic dissi- 
pation scale in the system, the scale of the disssipative 
stuctures is set by dispersion^ alone ( a Id collisionless soli- 
tons and shock waves). The dynamics of Alfven waves 
is sensitive to their sense of polarization, as well as to 
their propagation direction with respect to the ambient 
magnetic field. The principal results are: 

1. (quasi-), parallel initially linearly polarized waves 
evolve rapidly to form dissipative structures sep- 
arated by regions of linearly polarized Alfven 
waves. These structures, which are naturally quasi- 
stationary, extend the traditional genre of colli- 
sionless shocks to encompass the structures oc- 
curred via the competition of harmonic generation 
(i.e. nonlinear wave steepening) and collisionless 
(Landau) damping, rather than dispersion (i.e. as 
for solitons). Such structures are commonly oc- 
curring over a wide range of plasma (3, roughly: 
0.5 < (3 < 1.5. These structures are character- 
ized by fast phase changes (phase jumps). Namely, 
the magnetic field vector rotates first through some 
angle (typically, 50° — 90°), followed by rotation 
in the opposite direction through an angle larger 
than 180° to complete the total tt radian phase 
change. The phase diagram of the wave displays 
a distinct S-shape, which allows us to relate these 
dissipative structures to a class of S-type direc- 
tional/rotational discontinuities observed in the so- 
lar wind. Unlike conventional rotational disconti- 
nuities, for which = const, these KNLS rota- 
tional/directional discontinuities exhibit significant 
amplitude variation across the phase jumps. The 
amplitude of these variations depends on the an- 
gle of propagation and is largest for exactly parallel 
propagation. The KNLS rotational/directional dis- 
continuities have small velocity in the wave frame. 
Their directions are defined by the sign of the non- 
linearity coefficient Mi (the discontinuity propa- 
gates with the wave for AI2 > 0). The dissipa- 
tive character of these structures is revealed by 
their increasing width and decreasing amplitude. 
Dissipation is strong in this case, and energy de- 
cays roughly as E{t) oc l/(|M2|r). Since the 
KNLS directional/rotational discontinuities display 
a prominent S-shaped phase portrait, they can (po- 
tentially) be regarded as a signature of strong Lan- 
dau damping at work in satellite observations. 

2. (quasi-) parallel, initially circular polarizations do 
not form discontinuities. Since Landau damping 
preserves the initial symmetry of the harmonic 
spectrum, circular polarizations form intermittent, 
roughly arc-polarized waveforms which dissipate 
quickly, yielding purely circularly polarized waves 
at the lowest k. Such waves experience no nonlin- 
ear steepening and minimal damping. The differ- 
ent evolution of linear and circular polarizations is 
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due to specific coupling of the waves to collisionlcss 
damping, which is, in turn, sensitive to the (initial) 
symmetry of the harmonic spectrum. 

3. (quasi-) parallel, initially elliptically polarized 
waves represent an intermediate case between lin- 
ear and circular polarizations. Well localized, rota- 
tional discontinuities (similar to dissipative struc- 
tures) form. They are separated by large regions 
of linearly polarized waves. Amplitude variations 
across the discontinuities are small. These wave- 
forms are more weakly coupled to dissipation. The 
total energy of the waveforms decays as E{t) oc 
exp(-A|Af2|r), where A = - is the 
difference in energy between the two initially ex- 
cited harmonics. 

4. oblique, initially linearly polarized, longitudinal 
waves form KNLS rotational/directional disconti- 
nuities (dissipative structures) with opposite phase 
rotations, having different group velocities. Thus, 
these wave structures may interact to annihilate 
each other. As these waves are nonlinear ly coupled 
to the (large) ambient field B± = SgsinO, they 
are connected to dissipation more strongly than 
quasi-parallel waves are. Thus, the wave energy 
dissipate rapidly: E{t) oc exp(— |M2|t). Resid- 
ual, small-scale, weakly-nonlinear, weakly-damped 
Alfvenic perturbations may survive for a long time. 

5. oblique, initially circularly polarized waves and el- 
liptically polarized, longitudinal waves do not gener- 
ate dissipative structures. Despite this, they evolve 
to stationary waveforms with negligible amplitude 
and large phase variations. Their phase diagram 
displays a pure arc, with the arc-angle less than 
180°. Fast magnetic vector rotation along the arc 
at the discontinuity is followed by a part of an op- 
positely rotating circularly polarized wave. Thus, 
these structures have all requisites to be identi- 
fied with arc-polarized rotational discontinuities ob- 
served in the solar wind. This waveform experi- 
ences no damping and nonlinear steepening. How- 
ever, the energy content of this waveform is not 
small. 

6. oblique, transverse waves evolve to stationary arc- 
polarized rotational discontinuities with the arc- 
angle less than 180°. The energy content of this 
waveform is large. 

According to the discussion above, nonlinear Alfven 
waves subjcict to clamping on resonant particles c;volvc 
to narrow-spectrum waveforms with most energy con- 
centrated in lowest-A; harmonics. The initial spectrum 
symmetry (i.e. the wave polarization) is generally pre- 
served during the evolution. Based on the above, one can 
suggest that (given initial equal populations of isotropi- 
cally distributed circular and linear polarizations) quasi- 



parallel magnetic field fluctuations will consist of pre- 
dominantly circularly polarized waves and lower ampli- 
tude S-polarized KNLS directional discontinuities, while 
oblique perturbations are predominantly arc-polarized 
discontinuities, separated by pieces of oppositely circu- 
larly polarized waves. One should emphasize that all 
wave structures emerge naturally through the nonlinear 
evolution of simple, sinusoidally amplitude modulated 
Alfven waves with no a priori assumptions or special ini- 
tial conditions used. 

II. The dynamics of nonlinear Alfven wave trains 
subject to diffusive (collisional) damping differ signifi- 
cantly. There are two different regimes of weak and 
strong plasma coUisionality, controlled by the ratio of 
dissipation length L to the typical wavelength A, i.e. 
\L\ <C A and \L\ > A, respectively. The results are: 

1. The effect of collisions introduces nonlinear dis- 
persion into the dynamics, along with a nonlinear 
damping process. Thus it affects the direction of 
steepening, rate of spectrum modification, etc.. 

2. The sign of L defines the sign of nonlinear disper- 
sion which, in turn, is controlled by plasma /?, i.e. 
L (X {1 — (3). Thus, for /3 < 1 and /3 > 1, nonlinear 
steepening occurs at the rear and front edges of a 
train, respectively. 

3. The absolute value of L affects the dissipation rate. 
The damping rate has a maximum near |L| ~ A. In 
the weakly collisional regime, the dissipation pro- 
cess is slower on account of infrequent collisions. 
However in higher coUisionality plasma, the acous- 
tic mode is strongly damped, so it is no longer a 
plasma cigcnmodc (i.e. uj < 7. here 7 is the damp- 
ing rate). Thus, resonant energy transfer from 
Alfvenic to acoustic waves is suppressed. 

4. In a weakly collisional regime, a nonlinear wave 
evolves to a propagating shock wave (the direction 
of propagation coincides with the direction of wave 
steepening) and a rotational discontinuity or phase 
irregularities, depending on an initial polarization, 
i.e. whether the initial polarization is linear or cir- 
cular, respectively. The shock amplitude gradually 
decreases due to damping. The thickness of the 
shock is controlled by collisions, however the spa- 
tial size of the phase discontinuity/irregularities is 
larger and controlled by dispersion. In a strongly 
collisional regime, the wave evolution is much sim- 
ilar to the case of Landau damping. Obviously, the 
system lacks the characteristic dissipation scale in 
this limit, i.e. \L\ ^ A. 

5. The dissipation process is very sensitive (especially 
in weakly collisional plasma) to the value of /?, and 
sharply peaked at /3 equal to unity. If one assume 
P = P{R) is a function of the distance from the sun, 
then waves generated near the sun and propagating 
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outward will be strongly damped at the distance 
Rd defined by l3{Rd) = 1. One can, thus, expect 
increased plasma temperature in this region, as well 
as decreased magnetic field fluctuation level (i.e. 
characteristic amplitude of Alfven waves) beyond 
this point, i.e. at distances R > Rd- This process 
may be important, however, when the width of the 
dissipation region is large enough (comparable to 
or exceeding the typical wavelength). This is, in 
turn, controlled by the spatial profile of f3. 

III. Effects of finite Larmor radius impact the dy- 
namics of obliquely propagating waves, only. They are 
associated with small-scale gyro-averaging of fluctuating 
fields on a scale of Larmor orbits of plasma particles. 
These effects 

1. introduce nonlinear dispersion controlled by the 
plasma /?, Larmor scale pi, and angle of propaga- 
tion O. 

2. typically suppress the formation of arc-polarized 
discontinuities from initially circular polarizations, 
but do not change the qualitative evolution of lin- 
ear polarizations. 

3. affect the wave energy damping rate in high-/? 
plasma (/3 > 1) and may even change the sign of the 
dissipation in low-,/? plasma (/? < 1). The last case 
corresponds to an instability. Thus, as a result of 
gyro-motion of resonant particles interacting with 
a wave packet (in /3 < 1 plasma), an ion-cyclotron 
resonance occurs. It results in damping of low-A: 
harmonics and growth (instability) of high-A; har- 
monics, which interact with the tail of a particle 
distribution function. Thus, small-scale bursts may 
be generated in this regime. 

A more detailed investigation of obliquely propagating 
waves will be presented in future publications. 
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FIG. 1. Wave evolution from the DNLS with P — 0, Mi = 1, M2 — for linearly polarized, sinusoidal wave initial condition, 
(a) and (c) - wave profiles, (b) and (d) - jjhasc diagrams (hodographs) , (c) - harmonic spectra at r = 0,5,40. At r = only 
fc = 1, — 1 harmonics are excited. Note an irregular wave structure at late times (i.e. r = 40). 
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FIG. 2. Wave evolution from the DNLS with /3 = for circularly polarized, sinusoidally amplitude modulated initial condi- 
tion, (a) and (c) - wave profiles, (b) and (d) - phase diagrams (hodograplis), (c) - harmonic spectra at t = 0,8,40. At r = 
only k = —1, —2 harmonics are excited. Note the irregular wave structure at late times t = 40. 
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FIG. 3. Evolution for a parallel propagating, soliton initial condition for the KNLS with Landau damping for (3 = 1, Mi = 1.5, 
M2 = —1.63 at r = 0, 2, 25. (a) - wave packet spatial profile, (b) - hodographs, (c) - harmonic spectra, (d) - temporal evolution 
of wave energy, (e) - profiles of the wave magnetic fields bx and by at t = 25. 
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FIG. 4. Wave evolution from the KNLS with Landau damping for l3 — 1, AIi = 1.5, M2 = —1.63 for a parallel propagating, 
linearly polarized sinusoidal wave initial condition, (a) - wave profiles, (b) - hodograplis, (c) - harmonic spectra, (d) - temporal 
evolution of wave energy. At r = only fc = 1, — 1 harmonics are excited. Note the formation of localized, quasi-stationary 
waveforms {dissipative structures) with a characteristic S-shaped phase diagram and narrow harmonic spectrum. 
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FIG. 5. Detailed snapshot of a dissipative structure (or a KNLS rotational/directional discontinuity) at t = 15 in /3 = 1 
plasma, (a) - phase profile, (b) - wave amplitude profile, (c) - jjrofiles of the wave fields bx,by (note narrow spatial localization 
of the discontinuity) , (d) - hodograph (note a characteristic S-shape) , (e) -harmonic spectrum. 
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FIG. 6. Wave evolution for the KNLS with Landau damping for /3 = 1, Mi — 1.5, M2 = —1.63 for a parallel propagating, 
ciTculariy polarized sinusoidal wave initial condition, (a) - wave profiles, (b) - hodographs, (c) - harmonic spectra, (d) - temporal 
evolution of wave energy. At t = only k = —1, —2 harmonics are excited. Note the formation of an intermittent, arc-polarized 
discontinuity at r = 5, which disappears quickly due to Landau damping. Only k = —1 circularly polarized harmonic survives 
and experience no steepening or damping. 
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Landau, elliptical polarization 




FIG. 7. Wave evolution for the KNLS with Landau damping for /3 = 1, Afi — 1.5, M2 = —1.63 for a parallel propagating, 
elliptically polarized sinusoidal wave initial condition, (a) - wave profiles, (b) - hodographs, (c) - harmonic spectra, (d) - 
temporal evolution of wave energy. At t = only k = 1,-1 harmonics are excited and bi = 0.9, 6-1 = 1.1. Note the formation 
of rotational discontinuities (phase change equal to 7r radians) with weak wave amplitude variation across. Note also the smaller 
energy dissipation rate than for a linear polarization. 
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Landau, linear polarization, oblique propagation 




FIG. 8. Wave evolution for the KNLS with Landau damping for l3 — 1, Mi = 1.5, M2 = —1.63 for an obliquely propagating 
(O = 45°), /'tnear'/y polarized, longitudinal, sinusoidal wave initial condition, (a) and (c) - wave profiles, (b) and (d) - hodograplis, 
(e) - harmonic spectra, (f) - temporal evolution of wave energy. At t = only fe = 1,— 1 harmonics are excited. Note the 
formation of two KNLS directional discontinuities which interact and annihilate with each other (at r ~ 20) to yield a weakly 
dissipative, small amplitude residual wave structure. 
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FIG. 9. Wave evolution for the KNLS with Landau damping for /? = 1, Mi — 1.5, M2 = —1.63 for an obliquely propagating, 
circularly polarized, longitudinal, sinusoidal wave initial condition, (a) - wave profiles, (b) - hodographs, (c) - harmonic spectra, 
(d) - temporal evolution of wave energy. At t = only k = —1, —2 harmonics are excited. Note the formation of a stationary 
arc-polarized rotational discontinuity at approximately t = 25. 
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FIG. 10. Detailed snapshot of an arc-polarized rotational discontinuity at r = 40 for /3 = 1 plasma, (a) - phase profile, (b) - 
wave amplitude profile, (c) - profiles of the wave fields bx,by (note the narrow spatial locafization of the discontinuity and the 
negligible wave amplitude modulation), (d) - hodograph (note a characteristic arc-shape), (e) -harmonic spectrum. 
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Landau, linear polarization, oblique propagation, transverse wave 
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FIG. 11. Wave evolution for the KNLS with Landau damping for (3 = 1, Mi = 1.5, M2 = —1.63 for an obliquely propagating, 
trasversal, sinusoidal wave initial condition, (a) - wave profiles, (b) - hodographs, (c) - harmonic sjjcctra, (d) - temporal 
evolution of wave energy. At r = only k = 1,-1 harmonics are excited. Note the formation of a stationary arc-polarized 
rotational discontinuity at approximately r = 2. 
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Landau, harmonics 




FIG. 12. Temporal evolution of two initially excited harmonics for the KNLS for f3 = 1 for a (a) parallel propagating, 
linearly polarized wave (&_i and 61 harmonics overlaid), (b) parallel propagating, elliptically polarized wave, and (c) obliquely 
propagating (0 = 45°), linearly polarized wave. 
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FIG. 13. Evolution for a parallel propagating, soliton initial condition for the KNLS with coUisional damping for (5 = 0.99 
and K = 100 (k ~ ^typical corresponds to the strongest damping) at r = 0,4,25. (a) - wave packet spatial profile, (b) - 
hodographs, (c) - harmonic spectra, (d) - temporal evolution of wave energy, (e) - profiles of the wave magnetic fields bx and 
6„ at r = 25. 
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FIG. 14. Wave evolution for the KNLS with coUisional damping with /3 — 0.99 and k = 200 (weakly collisional regime) for 
parallel propagating, linearly polarized sinusoidal wave initial condition, (a) and (c) - wave profiles, (b) and (d) - hodographs, 
(e) - harmonic spectra. At t = only A: = 1, — 1 harmonics are excited. Note the formation of propagating shock waves (sharp 
fronts) and rotational discontinuities (wider wave- like forms). The noisy parasitic oscillations generated at earlier times (r = 1) 
are possibly due to modulation instability, which is stronger for smaller scale modulations. 
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CoUisional, linear polarization 




FIG. 15. Wave evolution for the KNLS with coUisional damping with /3 = 0.99 and k = 20 (intermediate collisional regime) 
for parallel propagating, /iraearij/ polarized sinusoidal wave initial condition, (a) and (c) - wave profiles, (b) and (d) - hodographs, 
(e) - harmonic spectra. At r = only k = 1,-1 harmonics are excited. The width of a shock wave is wider (as controlled 
by collisions), however the width of the rotational discontinuity is the same, i.e. it is controlled by the dispersion parameter, 
alone. 
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CoUisional, linear polarization 




FIG. 16. Wave evolution for the KNLS with coUisional damping with (3 = 0.99 and k, = 2 (strongly coUisional regime) for a 
parallel propagating, linearly polarized sinusoidal wave initial condition, (a) and (c) - wave profiles, (b) and (d) - hodographs, 
(e) - harmonic spectra. At t = only k = 1,-1 harmonics are excited. Note the remarkable similarity to the case of Landau 
damping. 
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FIG. 17. Detailed snapshot of a rotational discontinuity for an intermediate coUisional regime (k = 20) at r = 15 in /3 = 0.99 
plasma, (a) - phase profile, (b) - wave amplitude profile, (c) - profiles of the wave fields bx,by (note narrow spatial localization 
of the discontinuity) , (d) - hodograph (note a typical S-shape) , (e) -harmonic spectrum. 
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FIG. 18. Wave energy evolution: (a) - temporal evolution of wave energy im different coUisional regimes, (b) - wave energy 
content at times r = 5, 10, 20 as a function of plasma coUisionality, k. 
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CoUisional, linear polarization 




FIG. 19. Wave evolution for the KNLS with colhsional damping for /3 = 1.01 and k = —20 (intermediate coUisional regime) 
for a parallel propagating, linearly polarized sinusoidal wave initial condition, (a) and (c) - wave profiles, (b) and (d) - 
hodographs, (e) - harmonic spectra. At t = only k = 1,-1 harmonics are excited. Wavefront steepening occurs at the front 
face, thus it is controlled by the sign of k (the nonlinear dispersion effect). 
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FIG. 20. Wave evolution for the KNLS with conisional damping for /3 — 0.99 and k = 200 (weakly collisional regime) for a 
parallel propagating, circularly polarized sinusoidal wave initial condition, (a) - wave profiles, (b) - hodographs, (c) - harmonic 
spectra, (d) - temporal evolution of wave energy. At r = only k = —1, —2 harmonics are excited. A sharp front (shock) still 
forms, but phase irregularities form instead of a rotational discontinuity. 
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Finite Larmor, linear polarization 




FIG. 21. Wave evolution for the KNLS with Landau damping and gyro-kinetic effects for k?' < 0, |k| — 3000 (typical 
kpi ~ 0.002), and (3 = 0.9 for an obliquely propagating {O — 30°), linearly polarized sinusoidal wave initial condition, (a) - 
wave profiles, (b) - hodographs, (c) - harmonic spectra, (d) - temporal evolution of wave energy. At t = only k = 1,-1 
harmonics are excited. Observe that S-type rotational/directional discontinuities form. 
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FIG. 22. Wave evolution for the KNLS with Landau damping and gyro-kinetic effects for k?" < 0, |k| — 3000 (typical 
kpi ~ 0.002), and /3 = 0.9 for an obliquely propagating (0 = 30°), circularly polarized sinusoidal wave initial condition, (a) - 
wave profiles, (b) - hodographs, (c) - harmonic spectra, (d) - temporal evolution of wave energy. At r = only k = —1, —2 
harmonics are excited. Note that arc-type rotational discontinuities do not form. 
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FIG. 23. Wave evolution for the KNLS with Landau damping and gyro-kinotic effects for > 0, = 50 (large Larmor 
radius, typical kpi ~ 0.1), and l3 = 1.3 for an obliquely propagating (Q = 30°), linearly polarized sinusoidal wave initial 
condition, (a) - wave profiles, (b) - hodographs, (c) - harmonic spectra, (d) - temporal evolution of wave energy. At r = only 
A; = 1, — 1 harmonics are excited. 
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TABLE I. Summary of the numerical solutions. 
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